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Abstract. A series of sharp existence and uniqueness theorems are established for the 
multiple vortex solutions in the supersymmetric Chern-Simons-Higgs theory formalism 
of Aharony, Bergman, Jaferis, and Maldacena, for which the Higgs bosons and Dirac 
fermions lie in the bifundamental representation of the general gauge symmetry group 
U(N) x U(N). The governing equations are of the BPS type and derived by Kim, 
Kim, Kwon, and Nakajima in the mass-deformed framework labeled by a continuous 
parameter. 

1 Introduction 

It is well known that the presence of the Chern-Simons terms in field theories is essential in many 
areas of modern physics, especially for condensed matter systems [20- 22, 50J. More recently, the 
importance of the Chern-Simons terms in superstring theory and M-theory was explored in a general 
formalism by Schwarz |41j . in the context of some Lagrangian descriptions of superconformal gauge 
field theories which couple the Chern-Simons gauge and matter fields. This formalism was then 
made into fruition by Bagger and Lambert [HHZ] and Gustavsson |25| . which has since commonly 
been referred to as the Bagger-Lambert-Gustavsson (BLG) theory [9lll3|,[T8] . Shortly afterwards, 
Aharony, Bergman, Jafferis, and Maldacena (ABJM) developed an elegant bi-level Chern-Simons- 
Higgs theory [2], extending the BLG theory. In both the BLG and ABJM theories, topological 
solitons realized as the classical solutions of various BPS (named after the pioneering works of 
Bogomol'nyi [10] and Prasad-Sommerfield [40]) equations known as fuzzy funnels, domain walls, 
and vortices, are basic field configurations that describe M-branes. In the present study, we aim 
at establishing an existence and uniqueness theory for the multiple vortex solutions in a general 
class of BPS equations in the ABJM model. It is interesting to note that, in these superconformal 
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Lagrangian field theories [2"j[5HTll41"] , the dynamics of the gauge fields is exclusively of the Chern- 
Simons type, meaning that the usual quadratic kinetic terms involving the gauge field strength 
tensors, are missing [41], as in the earlier Chern-Simons-Higgs models of Hong-Kim-Pac [28] and 
Jackiw- Weinberg [30]. In these latter studies p8|30] and the subsequent development of the subject 
(cf. |16j and references therein), suitable six-order Higgs potential density functions have to be taken 
in order to achieve a BPS reduction. While, in such a situation, the BPS reductions can be made, 
the six-order potential terms lead to several complications or drawbacks of the models. The first 
one is that it is not clear whether topological solutions are uniquely determined by any prescribed 
distribution of vortices, as seen in the classical Abelian Higgs theory [31047]. although the existence 
of maximal solutions has been established [33]. The second one is the occurrence of non-topological 
solutions [29044] which are plagued by non-uniqueness and the issue of existence of solutions realizing 
arbitrarily prescribed vortices has only been tackled in some extreme cases [12[I14]. The third 
one concerns the solutions over a doubly-periodic domain resembling a vortex condensate/lattice 
structure [TJEIIIIE] and it is known that even in such a compact setting non-uniqueness occurs [35] 
and the interval ensuring the existence of a solution for the Chern-Simons coupling constant cannot 
be explicitly determined |114I45]. Furthermore, for non- Abelian BPS vortex equations, although 
the existence of topological solutions can be proved [51], the existence of vortex condensates has 
only been established for rank 2 gauge groups [39] and the existence of non-topological solutions 
is still unsettled [461153]. For the Chern-Simons vortices arising in the ABJM theory, however, 
we shall see through the present work that the situation is totally different. More precisely, we 
will develop a complete existence theory for multiple vortex solutions arising in the ABJM theory 
with the general gauge group U(N) x U(N) which gives us a unique topological solution for any 
prescribed distribution of vortices in the plane and a unique doubly periodic condensate solution 
under explicitly stated necessary and sufficient conditions involving several related physical coupling 
parameters. Furthermore, the model will be seen to be free of non-topological solutions. 

The governing BPS equations considered here for the multiple vortex solutions in the ABJM 
theory were discovered in the work of Kim-Kim-Kwon-Nakajima [33] who showed that for the 
bottom case when the gauge group is U(2) x U(2) the equations can be reduced to that of the 
Abelian Higgs equation [24,3U[38] which has been thoroughly understood [31U47] and that for the 
general case with the gauge group U(N) x U(N) (with N > 3) the equations are given as an (N— 1) 
by (N — 1) system of nonlinear elliptic equations of delicate structures. These equations arise in 
the so-called 'mass-deformed' theory labeled by a mass deformation parameter [i > 0. When fi = 
(mass deformation is absent), it has been shown in [33] that there is no finite-energy regular solution, 
and that finite-energy solutions can only be expected to appear when fx > (mass deformation is 
present). The aim of our work here is to establish an existence theory for such mass-deformed ABJM 
vortex equations. A similar structure has also been seen in the independent work of Auzzi and 
Kumar [4] for which some existence and uniqueness theorems have been obtained in [35] through 
exploiting appropriate variational methods. Enlightened by the success of that study, our technical 
strategy in the current study will again be to explore and unveil the underlying variational structures 
of the system. It is interesting to notice that, although the Higgs potential density function still 
contains six-order terms, which play a crucial role for generating non-topological vortices in the 
classical Chern-Simons-Higgs models |29|,I44] . such vortices are absent in the ABJM situation. 

The content of the rest of the paper is outlined as follows. In Section 2, we recall the multiple 
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vortex equations derived in the work of Kim-Kim-Kwon-Nakajima [33J in the ABJM model [2] 
and state our main existence results. The subsequent sections are then devoted to proofs of results. 
Specifically, in Section 3, we first apply a matrix decomposition procedure to unveil the variational 
structure of the system of governing elliptic equations. We next prove the existence of a solution 
by a somewhat indirect minimization approach, using the method of [31,37]. We then obtain the 
asymptotic decay properties of the solution on the full plane and calculate the anticipated quantized 
integrals which give rise to quantized vortex fluxes. In Section 4, we consider the compact case 
when solutions are doubly periodic. We shall mainly adapt the direct minimization method used 
in In Section 5, we present the limiting case a = which is of independent interest. We will see 
that such a limiting case allows us to state our necessary and sufficient conditions for the existence 
of a doubly periodic multiple solution explicitly. In Section 6, we reconsider the case when a > 
and present two concrete examples, with N = 3 and N = 4, respectively, as an illustration of the 
application of our results in the general situation (a > 0). We shall also describe how to use our 
results to estimate the dimension of the moduli space of the BPS equations in the ABJM model 
under consideration. 

2 Existence of vortices in the ABJM model 

Use fj,, v = 0,1,2 to denote the Lorentzian indices of the Minkowski spacetime M 2,1 of signature 
(- + +). Like the BLG model [5H71I25]. the ABJM model [2j[8] is formulated as a low-energy 
approximation of multiple M2-branes so that it is dual to M-theory on appropriate anti-de Sitter 
orbifolds. It is an N = 6 supersymmetric Chern-Simons theory with the gauge group U (N) x U (N) 
governed by the full Lagrangian density 

^ABJM = £ CS + £ kin ~ ^ferm ~ v ^ ( 2A ) 
in which the Chern-Simons Lagrangian £qq is given by 

C CS = ^ UXTr ( A »d v A x + ^A„A V A X - A„d v Ax - |i M i„i A l , (2.2) 

which describes two Lie algebra u(A r )-valued gauge fields, A^ and A^, with the Chern-Simons 
level (A;, —k); using 7^ to denote the Dirac matrices expressible in terms of the standard Pauli spin 
matrices, 7 = icr 2 ,^ 1 = cr 1 ,^/ 2 = cr 3 , the matter-kinetic Lagrangian >Cjjj n is given by 

Arin = " Tr (d^y\d»Y a ) + iTr (> A V Z^a) , (2.3) 

which couple four complex scalars (A = 1, 2, 3, 4) and four Dirac fermions V'A (A = 1> 2, 3, 4) in 
the bifundamental representations (N, N) of the gauge group so that the gauge-covariant derivatives 
take the form 

D^Y A = dpY A + \A il Y A — iY A A fl , 
D^Ya = d^Y A + \A^Y A -iY A A^ 

the Yukawa-like quartic-interaction potential density ^ erm is given by 

^ferm = ^Tr (y\y a ^ b -Y A Y^ B ^ + 2Y A Y^ A ^ 

-2Y\Y B ^ A H B - e ABCD Y^ B Y^ D + e ABCD Y A ^Y c ^) ; (2.5) 



.4 = 1,2,3,4; 



(2.4) 
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and Vo is the sextic scalar potential 
V = 



4tt 2 



— 5 Tc [6Y A YlY B Y\Y c Y^ - Y A Y\Y B Y^Y C Y^ - y\y a Y^Y b Y^Y c - 4Y a Y^Y c y\y b Y^ 



(2.6) 



With the Lagrangian density given by (|2.ip . the action 



^ABJM = J £ ABJM dx > ( 2 - 7 ) 
evaluated over the full spacetime is invariant under the N = 6 supersymmetry transformation [2)133] 



5Y A 

SlpA 

SA^ 
SA„ 



^ Ab iPb, 

--i^ABD^Y 13 + \ (-UJAB 

-1^ AB D^Y B + u BC (f3% c + <5jf 



Y c YqY b - Y b YqY c 



+ 2u bc Y b y\y c 



2tt 
" k 



IfiUAB + U) AB -i^lpAY ] B 



2tt 
k 



Y^uab + u^YX-iPb} , 



(2.1 



where ujab are supersymmetry transformation parameters satisfying 



,AB 



(uabY 



1 



ABCD , 



and /3^ B are given by the expressions 



pAB = ^LyiAy^yB] = *E (Y A Y ] C Y B - Y B YlY A 



(2.J 



(2.10) 



In the mass-deformed theory, we need to update the Lagrangian density (|2.ip by modifying the 
potential densities Vf erm and Vq following the recipe 



^ferm ^ ^ferm + AV ferm> V t->V + AV = V„, 



where 



AV ferm = Tr(^ A M^ B 
AV 



Tr | ^-Y A Y\Y B M%Yl - ^Y\Y A YlM B Y c + ^y\y a 



(2.11) 

(2.12) 
(2.13) 



in which /i > is a mass-deformation parameter which should not be confused with the Lorentzian 
index and M B = diag{l, 1, —1, —1}. We notice that, although the original Higgs potential density 
Vo given in (|2.6[) is purely sextic, the mass-deformed potential density V m obtained from adding 
the correction term AVo given in f)2. 13|) contains both quadratic and quartic terms, as that in the 
classical Chern-Simons-Higgs model [28l[30] . 

It is evident that the associated Euler-Lagrange equations of the mass-deformed action 

^mass-deformed ABJM = / ( £ CS + £ kin ~ ( F ferm + AF ferm) ~ V ™} dx ' ( 2 - 14 ) 
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are of course rather complicated. In the work of Kim, Kim, Kwon, and Nakajima [33J, it is shown 
that these equations in their static limit and in the absence of fermions may be reduced into the 
following remarkable BPS system of equations 



(d, \cr>„\v^ 


— n 




n t A =L 1 i — 1 0) 

— U, \Ji f= 1, l — X, 4), 


UQl i Lb y/J2 ' H' 1 ) 


— n 

— u > 


D Y 2 -is {Pf + ^Y 2 ) 


= o, 


D Q Y 3 - is/3i 13 


= o, 


D Y±-isPl* 


= o, 


ft 


= (3f=(3f+nY\ 


Pf 


= VY 3 , 




= HY\ 


Pf = P? = Pf = Pf 


= o, 


Pi C 


= o, (A^B^C^A) 



(2.15) 



where s = ±1 is a signature symbol to be specified later, coupled with the usual Gauss law 
constraints 



—B 

2vr 

^B 
2vr 



2tt 

K 

2n' 



'12 



J 







12 



-J 



(2.16) 
(2.17) 



where j° and j° are two associated matrix- valued conserved currents given by the expressions 

j, = i (y a D,y\ - D^Y A Y^ , (2.18) 
U = i MD,Y A - D,Y\Y A ) . (2.19) 



In [33], it is demonstrated that the energy of the mass-deformed ABJM model has the topological 
lower bound | 

(2.20) 



E> -n\Q + 2R 12 \, 



where the topological charges Q and R\ 2 are given by 

Q = Tr / j°dx, 



R12 = Tr J Jl 2 dx, 

for which the charge densities j° is defined in (|2. 18[) and J® 2 given by 

Jf 2 = i^DoY* - D Y x Yl) - \{Y 2 D Y% - D Y 2 Y}) 



(2.21) 
(2.22) 

(2.23) 



and the integration is carried over the full two-dimensional spatial domain, and the lower bound 
(|2.20p is attained by the solutions of the BPS system (|2.15[) coupled with the Gauss law constraints 
(|2.16p - ([2.17p so that s is determined by the condition 



\Q + 2R 12 \ = s(Q + 2R 12 ). 



(2.24) 



5 



Thus it is important to understand the solutions of (|2.15p - (|2.17p which will be our goal in the 
present work. 

To approach the system of equations (|2.15p - (|2.17p . Kim, Kim, Kwon, and Nakajima [33J take 
the following ansatz to represent the N x N matrices Y A = {Y^): 

Y i = ^$£fi> ^ = ^71h' y3 = ' y4 = ' (2 ' 25) 

where fi (i = 1, . . . , N — 1) are complex- valued functions and 

ai = yV + i - 1, i = l,...,N, (2.26) 

with a > a constant. Within this ansatz, the N x N matrix-valued 'magnetic' fields become 
diagonal whose entries are given by [33] 

= B^ = -2s^V + i - l){\fi\ 2 - |jU| 2 + 1), i,j = 1,...,N, (2.27) 

where the convention /o = Jn = is imposed. Using the complex-variable differentiation d = 
\{d\ — 182), d = \{d\ + 182), and eliminating the gauge fields from the equations, it is shown in [33] 
that the BPS system (|2.15p ~ (|2.17p is reduced into the following system of N — 1 (N > 3) coupled 
vortex equations: 

dd\n\h\ 2 = / i 2 ([2a 2 + l]|/ 1 | 2 -[a 2 + l]|/ 2 | 2 -l), (2.28) 

^lnl/,1 2 = ^(_^2 + ._ 1 ]| / ._ i |2 + J 2a 2 + 2 ._ 1 ]| / .|2_[ a 2 + .]| / . +i |2_ 1 ) ) 

i = 2, ...,N-2, (2.29) 
dBln\f N ^\ 2 = ^(-[ a 2 + N-2\\f N „2\ 2 +[2a 2 + 2N-3\\f N - 1 \ 2 -I) (2.30) 

away from the zero points of /1, . . . , /jv-i> which are known to be the vortex points of the system. 

Analyzing the structure of the system (I2.15P as in [31], it may be seen that the zeros of the 
fields fi, . . . , /jv_i are discrete and of integer multiplicities. Thus we can denote the sets of zeros 
of each fi by 

Z fi = {Pi,h---,Pi,ni}, i = l,...,N-l, (2.31) 

such that the number of repetitions of any point p among the set Zf { = {pi )S } (i = 1, ■ ■ ■ , N — 1) 
takes account of the multiplicities of the zero. 

We aim to prove that the prescribed sets of zeros given by (|2.3ip completely characterize the 
solution of (|2.28p - (|2.30p . To be precise we note that the problem may be considered either over 
the full plane M 2 under the natural boundary condition 

lim (591n|/,| 2 )(x) = 0, i = l,...,N-l, (2.32) 

\x\— Kx> 

or over a doubly periodic domain Q so that the field configurations are subject to the 't Hooft peri- 
odic boundary condition [26,48, 53j for which periodicity is achieved modulo gauge transformations. 
We shall establish that in both cases solutions exist and are unique. 
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Note that (pOg]) - ([2T3"0|) only make sense when N > 3. When N = 2 (the 'bottom' case with 
the gauge group U(2) x U(2)), the system is a single equation [33] : 

ddln|/| 2 = /, 2 ([2a 2 + l]|/| 2 -l), (2.33) 

which has been well studied and existence and uniqueness results have been obtained [3T1I471I48] . 
We now proceed to state our main results. 
Let R be the (N — 1) x (N — 1) tridiagonal matrix given by 



R = 

I 2a 2 + 1 -(a 2 + l) ... ... \ 

-(a 2 + l) 2a 2 + 3 -(a 2 + 2) ... 

-(a 2 + 2) 2a 2 + 5 -(a 2 + 3) ... 



'•■ -(a 2 + iV-3) 2a 2 + 2N-5 -(a 2 + N-2) 

\ ... -(a 2 + iV-2) 2a 2 + 2iV-3y 

(2.34) 

We will see later that the matrix R is positive definite. Denote the inverse of R by R~ 1 . We 
shall also see that all entries of R^ 1 are positive, i.e., (i? -1 )^- > 0, i, j = 1, . . . , N — 1. Write the 
eigenvalues of R as Ai, . . . , \n-i, and set Ao to be the positive quantity 

A = 2min{Ai,...,A JV _i}. (2.35) 

Our main results are collectively summarized as follows. 

Theorem 2.1 For any a > 0, /j, > 0, consider the system of multiple vortex equations f[2.28f l — f)2.30[) 
for the field configuration (/i, . . . , /jv-i) with the prescribed zero sets given by (|2.3ip such that each 
fi has ni arbitrarily distributed zeros pi t \, . . . ,Pi, ni , i = 1, . . . ,N — 1. 

(i) For the problem over the full plane M 2 , there exists a unique solution satisfying the boundary 
condition 

\fi\ 2 ^r u \x\^oo, i = l,...,N-l, (2.36) 
which realizes the boundary condition \2.32\) . where 

N-l 

r i =Y J (R' 1 h>0, i = l,...,JV-l. (2.37) 

3=1 

Moreover, this boundary condition is achieved exponentially fast at infinity, 

\\fi\ 2 ~n\ < C(e)e- 2 ^ £)VJr ^ as |s| -> oo, i = 1, . . . , N - 1, (2.38) 

where e G (0, 1) is arbitrarily small, C(s) is a positive constant depending on e, and Ao is defined 
by ([235D . 

(ii) Over a doubly periodic domain f2, the problem admits a unique solution if and only if 

N-l N-l 

nJ2( R ~ l h n 3<^\ n \Y,^ 1 ^' * = !,•••» ( 2 -39) 

3=1 3=1 
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hold simultaneously. 

(Hi) In both cases, there hold the quantized integrals 



f[^R ij \f j \ 2 -iUx = -^, 



i = l,...,N-l, 



(2.40) 



where the integration is evaluated either over the full plane M? or the doubly periodic domain Q. 

It can be checked that (|2.36|) and (|2.32p are equivalent. Thus, by virtue of (|2.37|) . we see that 
the so-called non-topological solutions [TSIHHIIHIIMIIII] do not appear in the ABJM model [2"II3"3"] 
considered here. 

We now proceed to compute the associated flux. 

As noticed in [33], the ansatz taken makes it consistent to assume that the gauge fields A\ and 
Ai (I = 1,2) are diagonal: 



Ai =Ai = diag{a i 1 ,...,of ,r }, 1 = 1,2. 
Thus, we may introduce the complex-valued variable 



(2.41) 



1 



A = -(Ax - L4 2 ) = diag <^ -(a\ - in*), -(of - i<) } = {h, . . . , b N } 



. , 1 j ' — ' ' . , i. ■ i ■ f . 1 1 i (2.42) 

so that the matrix-valued 'magnetic' field B = F\2 becomes 

B = F 12 = 9i^ 2 - d 2 Ax = -2i(dA - dA) = -2i diag{<%i - din, ...,db N - db N }, (2.43) 
which should not be confused with the group index used earlier. Hence, in view of (|2.15p . we have 

01n/i_i = (6i_i -h), i = 2,...,N, (2.44) 



away from the zeros of fi (i = 2, . . . ,N), where we have chosen s = — 1 for definiteness. From 
(I2.44p . we obtain 



i-1 



idd(]nf i - 1 + ]nf i _ 1 ) = (db i - 1 -db i „ 1 )-(db i -db i ), i = 2,...,N, (2.45) 



again away from the zeros of the functions. Comparing (j2.43j) and (I2.45|) . we see that B = 
diag{i?i, . . . , -Bat} over the same domain can be expressed as 

1 1 

5 i = ^_ 1 --Aln|/ 4 _ 1 | 2 = £ 1 --^A(ln|//), i = 2,...,N, (2.46) 



where B\ may be read off from (|237) . That is, B x = 2aV(|/i| 2 + 1). Consequently, we get 



N 



N-l 



Tr(B) = £ Bi = NBi - ~ £ (TV - i)A In |/,| 2 . (2.47) 

1=1 8=1 

In view of the equations (|2.28p ~ (|2.30p and the quantized integral formulas stated in (|2.40p . we 
see that over the doubly periodic domain Q the total 'magnetic' flux is 



/ MB) 

Jo. 



dx = 2Na 2 



N-l 



i + ^or-% 



8=1 



N-l 



N-l 



\n\ - ir {R'^im + 2tt (N - i)n u (2.48) 



i=l 



i=l 
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which is not quantized and depends on |fi|, unless a = 0. Here we have switched on the dependence 
on the signature symbol s = ±1 in the flux formula for generality. So the flux over the full plane 
M 2 diverges which leads to infinite energy as observed in [33]. Therefore, in order to avoid flux and 
energy divergence, it is of value to develop an existence theory for doubly periodic solutions when 
a > and of independent interest to spell out the existence theory separately when a = 0. Indeed, 
when a = 0, the matrix R simplifies itself considerably so that the results can be stated in concrete 
terms explicitly as follows. 

Theorem 2.2 For a = 0, \i > 0, consider the system of vortex equations (12.281) - (12. 30 j) for the field 
configuration (/i, . . . , /jv-i) with the prescribed zero sets given by (|2.3ip such that each fi has n« 
arbitrarily distributed zeros p^i, . . . ,Pi )Ui , i = 1, . . . , N — 1. 

(i) For the problem over the full plane M 2 , there exists a unique solution satisfying the boundary 
condition 

\fi\ 2 ^(N-i), \x\^oo, i = l,...,N-l. (2.49) 
Moreover, this boundary condition is achieved exponentially fast at infinity 

\\fi\ 2 - {N-i)\ < C(e)e^ (1 - £) ^° N , i = l,...,N-l, (2.50) 

where e G (0, 1) is arbitrarily small, C{e) is a positive constant depending on e, and Ao is defined 
by (^351) . 

(ii) For the problem over a doubly periodic domain Q,, there exists a unique solution if and only 
if the conditions 

N-l j 

-E n ' <» 2 \n\(N-i), i = l,...,N-l, (2.51) 
j=i 3 i=i 

hold simultaneously. 

(Hi) In both cases, there hold the quantized integrals 

(|/i| 2 -|/ 2 | 2 -l)dx = -^, (2.52) 

J (-[^-l]|/ i _ 1 | 2 + [2i-l]|/,| 2 -i|/ m | 2 -l)dx = -^, i = 2,...,N-2, (2.53) 

j (-[N-2]\f N ^\ 2 + [2N-3}\f N ^\ 2 -l)dx = -^l, (2.54) 

where the integration is evaluated either over the full plane M 2 or the doubly periodic domain fL 

To see the problem more transparently, we reformulate the system of equations (|2.28|) f)2.30[) 
and Theorems I2.1H2.2I in terms of a new family of parameters and variables. For this purpose, we 
denote N — 1 = m, A = 4/i 2 , and set 

iii = ln|/j| 2 , i = l,...,m. (2.55) 
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Then, in view of the equations (|2.28p - (|2.30p and the zero sets (|2.3ip we see that u\, . . . , u m satisfy 
the equations 



A Ul = A ( [2a 2 + l]e Ul - [a 2 + l]e" 2 - l) + 4tt ^ 5 Pls 



(2.56) 



Am = A (-[a 2 + i-l]e u *~ 1 + [2a 2 + 2i-l]e u ' -[a 2 + i]e u * +1 - l) +4ttJ^5 Wi ., (2.57) 

8=1 

i = 2, . . . m — 1, 

Au m = \{-[a 2 + m-l]e Um ~ 1 + [2a 2 + 2m - l]e" m - l) + 4vr ^ <y pmiS . (2.58) 



Let 



u 
1 

U 



(Ul, ■ ■ ■ ,Um) T , 

(i,...,ir, 

(e u \...,e" m ) T 



5 = 1 8=1 



(2.59) 



Then the equations f|2.56j) — ()2.58j) can be recast into its vector form 

Au = A (RU - 1) + 4vrs. 



(2.60) 



Theorems I2.1H2.2I will be established through proving the following results for the equations 



Theorem 2.3 For any a > 0, A > 0, consider the system of equations ()2.56j) ()2.58[) . 

(i) There exists a unique solution overM? satisfying the boundary conditions 



Ui 



In [ ^(R- 1 ) 



= lnrj, |x| — > oo, i = l,...,m. (2.61) 
Moreover, this solution satisfies the following exponential decay estimate at infinity: 

m 

Y,( u i(x) - hir,) 2 < C(e)e-^- £)VJ ^^, (2.62) 



i=l 



where e 6 (0, 1) is arbitrarily small, C{e) is a positive constant depending on e, and Ao is defined 
by flZMD - 

(ii) For the problem over a doubly periodic domain £1, a solution exists if and only if the following 
m conditions 



4vr^(i? x ) iS nj < \\tt\J2(R i = l,...,m, 

j=i i=i 

hold simultaneously. Besides, if a solution exists, it must be unique. 



(2.63) 
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(Hi) In both full plane and periodic domain cases, there hold the quantized integrals 

J R ^ Uj - 1 j dx = -^T' » = 1, ■ ■ ■ , ™, (2.64) 

evaluated over the corresponding domain of consideration. 
When a = 0, we have the following explicit results. 

Theorem 2.4 For a = 0, A > 0, consider the system of equations (I2.56H ( 12.581) . 

(i) There exists a unique solution overW? satisfying the boundary conditions 

Ui — >• ln(m — i + 1), \x\ — y oo, i = l,...,m. (2.65) 
Moreover, this solution satisfies the following exponential decay estimate at infinity: 

m 

(ui(x) - ln(m - i + l)) 2 < C(e)e- ( - 1 - £)VJX ° lxl , (2.66) 

i=i 

where e E (0, 1) is small, C(e) > depends on e, and Ao is defined by (|2.35[) . 

(ii) There exists a solution over a doubly periodic domain fi if and only if the conditions 

m 3 

4vrJ^-^n/ < X\n\(m-i + l), i = l,...,m, (2.67) 
j=i 3 i=i 

are fulfilled simultaneously. Besides, if a solution exists, it must be unique. 
(Hi) In both cases, the quantized integrals 

j (e Ul -e U2 -l)dx = (2.68) 

/4-7T77 ' 
(-[%- lje^" 1 + [2i- l]e u > - ie u ^ - 1) dx = -±, i = I, . . . ,m - 1, (2.69) 

J {-[m- l]e u — 1 + [2m - l]e Um - 1) dx = (2.70) 
are uaZid ower the corresponding domain of the problem. 
In the subsequent sections, we prove Theorems I2.3H2.41 

3 Variational principle and solution to planar case 

In this section we establish the existence and uniqueness results for a solution of (|2,56p - (|2.58p over 
the full plane and derive the stated decay estimates for the solution. Unlike the problems studied 
in [35H37] which can be readily formulated variationally, the problem here needs more elaboration in 
order that its hidden variational structure be unveiled. For this purpose, we shall rely on the well- 
known Cholesky decomposition theorem for positive-definite matrices. The variational structure to 
be recognized will allow us to prove the existence of a solution over the full plane as well as over a 
doubly periodic domain, although the present section is devoted to the planar case. Below we split 
our study into a few subsections. 
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3.1 Cholesky decomposition for the matrix R 

First, we observe that the symmetric matrix R is positive definite. In fact, for any a > 0, it is easy 
to check that each leading principal minor of R is positive. That is, 



Ri = det([Rij] M ) > 0, i,j = l,...,l, l = l,...,m. (3.1) 

By the Cholesky decomposition theorem [23] the matrix R can be expressed as the product of a 
lower triangular matrix L and its transpose, R = LL T , L = (Lij) mxm . Indeed, using the iteration 
scheme presented in [23], that is, 



La 



L 



& i = 2,.. .,m, 

I —i 
jRii - Y, L 1p i = 2,...,m, 

3-1 

T 1 ) i=j + l,...,m,j = 2,...,m, 



(3.2) 



we have 



Li 



Lu—i 
Lij 



\/2a 2 + 1 = VHi, ^21 
.( a 2 + i _!) 



(" 2 +l) 



Ri-i ' L " 
0, l<i<t-l, i = 2,. 



Rj 

Ri-i' 
. , m. 



i = 2 



, ... ,m, 



1 



Here and in the sequel, we follow the convention i?Q = 1. We have 



/ V 7 ^ 

q 2 +l 









(a^ + 2) A /f 








V 



-(a 2 + m- T 



Rm-2 
Rm-1 



o \ 







Rm-1 / 



(3.3) 



(3.4) 



Furthermore, a simple calculation enables us to find the lower triangular matrix L 1 with 



(a 2 +j)---(q 2 +i-l)R 

Ri-lRi 
R i~l j — 2 

Ri ' ' ~ z ' 



■i-i 



1 < j < i — 1, i = 2, . . . ,m, 



, m. 



(3.5) 



By the expression of L 1 we can compute the inverse of f? by the formula 

R- 1 = {LL T )- 1 = (L^YL- 1 , 



(3.6) 



from which we can see directly that, for any a > 0, all entries of R 1 are positive: (R > 
0, i, j = 1, . . . , m, as claimed earlier. 
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3.2 Variational formulation 

Following [51] . we introduce the background functions 



^ln(l + v\x - p it s\ 2 ), i = l,...,m, 



(3.7) 



s=l 



with ^ > being a parameter which should not be confused with the Lorentzian index used earlier. 
We see that u l satisfy 



4;y 



s=l 



•[{v+\x-pi 



\2\2 



i = 1, . . . , m. 



(3i 



It is easy to see that 



gi G L(I 2 )nL 2 (I 2 ) 



/ 



jdx = 47rnj, % = 1, . . . ,m. 



(3.9) 



Since for any a > each entry of i? is positive, we have 



ri = ^2(R-% >0, i = l,... 

3=1 



, m. 



(3.10) 



Set r = (n, . . . ,r m ) T and Uj = + lnr* + «», i = 1, . . . m, v = . . . , v m ) T , g = (g%, . . . ,g m ) T , 
U = (e n o+^i, . . . ,e u o+ v ™) T . Then the equations <[2T56|> - ff2T58|l over M 2 become 



A Vl = A ([2a 2 + l]ne u o +vl - [a 2 + l]r 2 e< +V2 - lj + 9l , (3.11) 
Avi = \(-[a 2 + i-i\r i -ie< 1+Vi ~ 1 + [2a 2 + 2i - l]r t e u ° +v * 

-[a 2 +i]r i+1 e< 1+Vi + 1 - l) + 9i , * = 2, . . . m - 1, (3.12) 
Aw m = A (-[a 2 + m - l\r m - X e<~ 1 +v ™~ 1 + [2a 2 + 2m - l]r m e u ™ +,,m - l) + g m . (3.13) 
Or in its equivalent vector form, we have 

Av = Ai?diag{n, . . . , r m } (U - 1) + g, (3.14) 
since R diagjri, . . . , r m }l = 1 by the definition of r\, . . . , r m . Now we use the notation 

w = (u>i, . . . , w m ) T , h=-L _1 g, h = (hi, . . . ,h m ) T (3.15) 
and the transformation 

/ (L- l ) n vi \ 



w = L v 



\(L l ) m \Vl H h (L 1 )mmV m J 



(3.16) 
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or 



v = Lw 



/ Luwi \ 

L21W1 + L22W2 



(3.17) 



y-^'mm-l'i'm-l ~i~ ^ram^m J 

Then the equations f|3. 1 1 P — (|3. 13j) take the form: 



Aw\ = A \ L\\r\ 
Awi = \\LaTi 



+ L>2\T2 



JIQ+L21W1+L22W2 _ -y 



+ hi 



(3.18) 



+Li+uri + i 



Aw m = A {L mm r m [ e < n +L mm -iw m . 1+LmmWm _q + hm) _ 
Or in its vector form, we have 

Aw = A (£ T diag{n, • • • , r m } [U - 1] + h) . 



+ hi), i = 2,...,m-l, (3.19) 

(3.20) 



(3.21) 



It can now be checked to see that the equations (|3.18p - (|3.20p or (|3.2ip are the Euler-Lagrange 
equations of the functional 

J(w) = — V / \Vwi\ 2 dx + y^ h iWl dx+ n (e u o +LllWl - e u o - LnwA dx 



i=2 



+ J2 ri Lui+Lu-im-i+LuVi _ e w _ l Lii _ xWi -i + Luwi]) dx 



(3.22) 



This is the variational principle we have aimed to unveil. 

To facilitate the computation and analysis, it will be technically more convenient to rewrite the 
functional ()3.22j) as 

J(w) = i- V / \Vwi\ 2 dx + [ ne u ° (e LnWl - 1 - L n wA dx 

+ ^ / r ie u o ^Lu-iw^+LuWi _ 1 _ [L^iti;^! + dx 



i=2 
m—1 



e"o - 1 



i+li' i+1 



+ m L ' 

[e u o - 1] + /i m ) w m da;, 

Jk 2 



' i + l 

e u o -1 



+ hi Widx 



(3.23) 



which allows us to approach the problem in a similar manner as in |3iy47j for the scalar situation, 
as we will do in the following. 
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To proceed, we can compute to get 
(D/(w))(w) 



I+L22W2 



ll ) widx 



= I J] f \V Wi \ 2 dx + f finne^ [e illWl - l] + L 21 r 2 e< [e L ^~ 
A ~[ Jr. 2 Jr 2 v 

m— 1 . 

+ V / (Lun^o (eLu-im-i+Luwi _q + Li+un+ie ^ +1 r e L i+u w i+ L i+li+1 w i+1 _ -n\ w . dx 

~o JR 2 ' 

{L mm r m e u ° [e L ^-iw m -i+L mm w m _ x ]) Wmdx 

+ E L ( L < 



m—1 



e u ° - 1 



i+l 

e"o - 1 



+ hi Widx 



+ / {Lmmr m [e U ™ - l] + /i m ) W m dx 
JR 2 

= |V / |V^| 2 dx+ / n f e "o +I/llWl - 1 + fr^ Lnioidx 

A ~[ Jr 2 Jr 2 v y 

m „ 

+ Y n Lui+L^w^+L^ _ 1 + f li ) (Lu^Wi^ + LuWi) dx, 

where /ij's are some linear combinations of h^s, also of g^s. More precisely, 

hi = Y~.^2( Rl )ji9j, i = l,...,m. 



(3.24) 



(3.25) 



Using the invertibility of the transformations (I3.16|) and (|3.17p . we see that there exist some 
positive constants c\ and C2 such that 

m m m m m m 

ci Y, v l - XX - ° 2 J2 v l c iY,\ Vvi \ 2 -Y,\ Vwi \ 2 - c2 J2\ Vvi \ 2 - ( 3 - 26 ) 

i=l i=l i=l i=l i=l i=l 

Therefore, from (|3.17p . (|3.24p and (|3.26p . we can obtain 

m . m „ 

(£>I(w))(w) >C J2 |V^| 2 dx + ^ / r, (e M o+^ - 1 + ^dx, (3.27) 
i=i ^ R2 i=i ^ R2 

where Co is a positive constant. 

To deal with the second term on the right hand side of (|3.27p . we follow the approach of |31| . 
We just need to estimate a typical term of the following form 

M(v) = [ (e Uo+v - 1 + h) vdx. (3.28) 



It is easy to see that M(v) can be decomposed as M(v) = M(v + ) + M(-v-) where v + 
max{i>, 0},v_ = max{— v, 0}. 

From the elementary inequality e* — 1 > t, t G M and the fact no, h € L 2 (R 2 ), we have 

1 



M(v+) > I (u + v + + h)v + dx > - I v+dx-Ci, 



(3.29) 



15 



where and in the sequel we use C to denote a generic positive constant. 

t 

l+t 



To estimate M(-v-), we note the inequality 1 — e * > y4r, t > 0, and obtain 



M(-v. 



> / 

JR- 



1 - e u °~ v - - h ) v-dx 
1 - h - e"° + e uo [1 - e' 
l-h-e u °+e u ° 



V-dx 



1 + V- 



V-dx 



[1 + v 
1 - h 



l-h-e U0 

2 



+ e Uo v 



1 + v_ 



-dx 



v 



> 



> 



1 + v. 



1 + V- 



-dx + 



dx + 



1 -e 



1 + V_ 



-dx 



1-e* 



1 + v. 



-dx 



+ 



r dx - C 2 , 



(3.30) 



where we have used h < |, assured by taking i/ sufficiently large, and the fact that both e u ° — 1 
and h belong to L 2 (M. 2 ). 

Then we see from (|3.29p and (|3.30j) that 



M(v) > - [ 



rdx - C. 



(i + M) 2 

To proceed further, we need the following standard interpolation inequality over W 1,2 (M. 2 ): 



(3.31) 



/ t> 4 dx<2/ v 2 dx [ |Vw| 2 dx, V-u G W 
Jr 2 Jr 2 Jr 2 



(3.32) 



Using (|3.32|) . we have 



\v \ 2 dx 



v 



l + \v 

|2 



< 

< 4 
1 

< - 



v 



(i + H) 



(1 + |u|)|u|dx 

^■dx / (\v\ + b| 2 ) 2 dx 
' Jr 2 



(1 + M) 2 



dx / |w| 2 dx ( / |Vv| 2 dx + 1 ) 
Jr 2 \Jr 2 J 



v\ 2 dx) +C 



v 



(i + M) 2 



dx 



+ 



/ |Vv| 2 dx 
Jr 2 



+ 1 



(3.33) 



which implies 



\vh < C 



(i + \v\y 



; dx+ / |V?;| 2 dx + 1 ] 
Jr 2 J 



(3.34) 
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where and in the sequel we use || • [L to denote the norm of the space L' P (M?). 
From (ET2T1) and (pOT]) . we obtain 

(D /(w))(w) > c 2 g I (|v„/ + j d, - c 3 . 

Then it follows from (13^61) . (13321) and (^351) that 



(3.35) 



(£>J(w))(w) > C 4 ^ 11^11^/1,2^2) -C 5 . (3.36) 
j'=i 

By the coercive lower bound in (|3.36p . we can show that the functional / defined by (|3.22p 
admits a critical point. In fact, by (|3.36p . we can choose £ > such that 

mi{DI(w)(w) | 1^11^/1,2^2) = £} > 1. (3.37) 

Since the functional / is weakly lower semi-continuous on W 1,2 (M. 2 ), the minimization problem 

r/o = inf {/(w)| 1^11^/1,2^2) < £} (3.38) 

admits a solution, say, w. Then, we can show that it must be an interior point. We argue by 
contradiction. Assume that ||w||wi,2(k2\ = £. Then 



a8 /((l- t )w)-/(w) = d f((1 _ wu = < _! (3 . 39) 

Therefore, if t > is sufficiently small, letting w — (1 — t)w, we can obtain 

I(w f ) < I(w) = r/o, 11^*11^1,2^2) = (1 - t)£ < £, (3.40) 



which contradicts the definition of r/o. Hence, w must be an interior critical point for the problem 
(|3,38p . As a result, it is a critical point of the functional /. Since the functional / is strictly convex, 
this critical point must be unique. 

3.3 Asymptotic behavior of the solution at infinity 

Here we study the asymptotic behavior of the solution obtained above. Noting that w S iy 1,2 (lR 2 ), 
by the well-known inequality 

\\e v - l\\l < C x exp(C 2 \\v\\ 2 wl , HR2) ), Vv G W 1,2 (R 2 ), (3.41) 

where Ci,C% are some positive constants, we see that the right-hand sides of the equations (]3.18p ~ 
(|3.20|) all belong to L 2 (R 2 ). Using the standard elliptic L 2 -estimates, we have Wj £ W 2 ' 2 (M?), 
which implies Wj — > as \x\ — > oo, j = 1, . . . , m. By the transformation (|3.17p . we see that Vj — > 
as |x| —> oo, which implies the desired boundary condition Uj — > lnrj as \x\ —> oo, j = 1, . . . ,m. 

Next we show that | ViUj| — > as \x\ — > oo, j = 1, . . . , m. A typical term of the right hand sides 
of (|3.18p - (|3.20p can be rewritten as 



e u o+ L n-i w i-i+ L 33 w 3 _ I 



^ e U — ij e L jj-l W j-l+ L n W j _|_ (,.'•;; I "'I : • _ \j ; (3.42) 
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which belongs to L P (M 2 ) for any p > 2 due to the embedding W 1,2 (R 2 ) C L P (R 2 ) and the definition 
of Uq. Therefore all the right-hand-side terms of (|3.18p - (|3.20p lie in L P (M. 2 ), for any p > 2. Then 
the elliptic L p -estimates imply Wj G W 2,P (M. 2 ) for any p > 2, j = 1, . . . , m. Consequently, we have 
\Vwj\ — > as \x\ — > oo, j = 1, . . . m. That is, \Vuj\ — > as \x\ — > oo, j = 1, . . . m. 

Now we establish the exponential decay rate of the solutions at infinity. To this end, we consider 
the equations (|2.56p - (|2,58p or (|2.60p over an exterior domain 

D p = {x G M 2 | |x| > p) , (3.43) 

where p > satisfies 

p > max {\pi )S \ | i = 1, . . . ,m, s = 1, . . . ,m}. (3.44) 

For convenience, we consider the system of equations (I2.60P over D p . Recall that r = 
Hence we may rewrite (|2.60p in D p as 

Au = \R(U - r) = Ai?v + \R(U - r - v) , (3.45) 

where vector v is to be determined. 

Since the matrix R is positive definite, there is an orthogonal matrix O such that 

O t RO = diag{Ai,...,A m }, mm{Aa,...,A m }>0. (3.46) 

Now apply T in (|3.45p and set 

u = T v, v = (ui - lnri, . . . ,u m - lnr m ) T . (3.47) 

Then we have 

ACi = Adiag{Ai, . . . , A m }Ci + XO T R(XJ - r - v). (3.48) 

Note that, since U — > r as \x\ — > oo, we have U — r = E(x)v where E(x) is an m x m diagonal 
matrix so that E(x) — > I m (the m x m identity matrix) as \x\ — > oo. This observation leads us to 
rewrite (|3.48|) as 

ACi = Adiag{Ai, . . . , A m }£i + XP(x)u, (3.49) 
where P(x) is an m x m matrix which vanishes at infinity. As a consequence of (|3.49p . we obtain 

A|u| 2 > 2£i r Au > AA |u| 2 - b(x)\u\ 2 (3.50) 

with b(x) — > as |x| — > oo. 

Then, for any e G (0, 1), we can find a suitably large p e > p such that 

A|u| 2 > (l - |) AA |u| 2 , x£D P£ . (3.51) 

We will use a comparison function, say r/, of the form 

77 = Ce~ a W, \x\ > 0, C, a G R, C, a > 0. (3.52) 

Then Ar/ = a 2 r\ — ^V- Thus, in view of (|3,5ip . we have 

A(\u\ 2 - V ) > AA |u| 2 - ex 2 //, \x\>p £ . (3.53) 
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We take the obvious choice a 2 = (l — §) AAo which gives us A(|u| 2 — rj) > a 2 (\u\ 2 — 77), |x| > p £ . 
Choose C in (|3.52p large so that |u| 2 — rj < for |x| = p £ . Hence, using the fact that |u| — > as 
\x\ — > 00 and the maximum principle, we see that |u| 2 < rj for \x\ > p £ . So the estimate 

|£i| 2 < C&e-^^W, \x\ > p £ , (3.54) 

follows since (1 — e) 2 < (l — | ) for any e 6 (0, 1). Therefore the desired exponential decay rate 
(f2T62|) is established. 

To get the quantized integrals, we need to establish the exponential decay rate for the derivatives. 
Let d denote any of the two derivatives d\ and 82- Define 

v = (du 1: . . . ,du m ) T , P = diag{ri,...,r m }, Q = diag{e ui - r h . . . , e Um - r m }. (3.55) 

Then differentiating ()2.60p in D p , we have 

Av = XRPv + XRQv. (3.56) 

Let O be as before and set 

Fv = Ou, f = v T Pv. (3.57) 
Then by f|3.56j) and the fact that r\, . . . ,r m > 0, we obtain 

A/ > 2v T PAv = 2\v T PRPv + 2Xv T PRQv 

> AA |u| 2 + 2Av T RRQv = AA v T P 2 v + 2\v T PRQv 

> XX o r f-b(x)f, (3.58) 

where ro = min{ri, . . . , r m } and b(x) — > as \x\ — > 00. Then, as discussed previously, we can 
conclude that for any e G (0, 1), there is a positive constant C(e) > 0, such that 



/ < C(e)e- (1 - e)vTO ° |a!| , (3.59) 
when \x\ is sufficiently large. Thus we get the following exponential decay rate near infinity: 



l v ^0*0l 2 < C(e)e-( 1 - £ )^^°l a: l. (3.60) 



i=i 



We can now calculate the quantized integrals (|2.64p stated in Theorem 12.31 for the planar case. 
Using (|3,7p . ()3.8p . and the exponential decay property of | Vita's in (|3.60p . we conclude that 
|Vfj|'s vanish at infinity at least at the rate |x| -3 . Thus, using the divergence theorem, we have 

Avi dx = 0, i = l,...,m. (3.61) 

Consequently, integrating the equations (|3.1ip - (|3,13p over M 2 , and applying (|3.9p and (|3.6ip . we 
obtain the desired results stated in (|2.64p . 

We next turn our attention to the compact case. 
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4 Doubly periodic case 

In this section we consider the equations (|2.56j) — (|2.58|) over a doubly periodic domain Q. 
Let Uq be a solution of the problem (see [3]) 

A4 = -^i + 47T^<5 p . s; xefi, i = l,...,m. (4.1) 

Set in = u l + Vi,i = I, . . . ,m. Then equations (|2.56p - (|2.58p become 

A Vl = \([2a 2 + l]e< +v ' -[a 2 + l]e< +v * -l) (4.2) 



Avi = A (- [a 2 + i - 1] e n o _1 +^-i + [2a 2 + 2* - l] e u o +Vi - [a 2 + i] e "o +1 +^+i - l 

+ 7^f' i = 2 '--- m - 1 ' ( 4 - 3 ) 

At; m = A(- [a 2 + m-l]e<~ 1+v ™-i + [2a 2 + 2m - l] e u ™+^ - l) + (4.4) 
or equivalently in its vector form 

Av = \(RU - 1) + ^n, (4.5) 

withU= (e u o+^,..., e ^+^) T ,v = (vi,...v m ) T ,l = (1, . . . , 1) T , n = (n u . . . ,n m ) T . 

Integrating the equations (|4.2p - (|4.4p or (|4.5p over f2, we can obtain the natural constraints 

R / Udx = - — n. (4.6) 
Jn A 



We may rewrite (|4.6p more conveniently as 



/ 



4tt 

Udx = \n\R~ 1 ! - — R- 1 !! = K, (4.7) 
n A 



or in its component form 

„ m . m 

/ e< +Vi dx = V(fr% - V(tt : ')y n i = » = 1, ■ ■ • , m, (4.8) 

where and in what follows we use the notation K = (K\, . . . ,K m ) T . Therefore, we see that if a 
solution exists, then ifj > 0, i = 1, . . . , m. As a result, we get the necessity of the condition (|2.6ip . 

In what follows we show that the condition (|2.61|) is also sufficient for the existence of a solution 
to (j4.2p - (j4.4p by variational methods. 

We will work on the Sobolev spaces W l,2 {£i), which is composed of scalar- or vector- valued 
O-periodic L 2 -functions whose derivatives also belong to L 2 (Cl). For the scalar case we have the 
decomposition 

W 1 > 2 (Sl)=R@W 1 ' 2 (n), (4.9) 



where 



w^ 2 (n) = {we w l ' 2 {Q) 



wdx = } , (4.10) 
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is a closed subspace of W 1,2 (Q). Then for any u £ W 1,2 (£l), we have 

u = u + u, u£R, iteW 1 ' 2 ^). (4.11) 
We will use the well-known Trudinger-Moser inequality [3"lll9| 

J e u dx < Cexp (r^ J |Vn| 2 dx^, VuG ^(fi), (4.12) 

which is important for our estimate, although the analysis does not depend on the optimality of 
the embedding constant. 

As in the planar case, to formulate the problem in a variational structure, we use the transfor- 
mation (|3.16|) or (|3.17p . Then the equations ()4.2p - (j4.4p become 

A Wl = X (L n e u o +L ^ m + L 21 e u ° +L2lWl+L22W2 - bA , (4.13) 

A Wl = X ^L iie u o +L "' im - 1+Lum + L l+li e u ° +1+Ll + um+Ll + u + im+1 - b^j , (4.14) 
i = 2, . . . , m — 1 

Aw m = X (L mm e u ™ +Lmm - lWm - 1+LmmWm - b m ) , (4.15) 
whose vector form is 

Aw = A(L r U-b), (4.16) 

where 

b = (h, . . . , b m f = L-H - jr^L^n, (4.17) 

As before, we can check to see that the equations (|4.13p - (|4.15p are the Euler-Lagrange equations 
of the functional 

J(w) = — V / |V^| 2 dx+ / e^' + Ve^ 1 ^'* di-V / hwidx. (4.18) 
2A . =1 Jn J n \ i=2 J i=1 Jn 



We shall now engage ourselves in a direct minimization procedure, initiated in [35], to find a 
critical point of the functional (|4.18p . 

When w E W 1,2 (£l), by the Trudinger-Moser inequality (|4.12p . we see that the functional 
defined by (|4.18p is a C 1 -functional and lower semi-continuous with respect to the weak topology 
of W 1 ' 2 ^). 

Using the decomposition formula (|4.1ip . we obtain 



-. m „ / m \ 

j( w ) V || Vm \\l = / I e u ° +LllWl + V e ^h+Lii-im-i+L»Wi \ 

2X tl " J n\ J 

m „ 

= V / e"o +1,l dx-K T v 
i=i Jn 



dx — |Sl|b T w 



rn ,. m 

/ e^^dx-^KiVi, (4.19) 
i=i ^ n i=i 
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where we have used |f2|b r w = K r v in view of (|4.7p and (|4,17p . 
By Jensen's inequality, we see that 

J e «S+**+s* dx > \n\ exp ( y (u* + ^ + ^) dx^) 

= |«| exp ( y 4dx^ e 8 * = trie 2 *, i = l,...,m. (4.20) 

Using the condition (|2,63p . we have K{ > 0, i = 1, . . . , m. Then, combining (|4,19p and (|4.20p . 
we have 



1 m „ m 
i=l Jn i=l 



i=i 



Hence, from (|4.2ip . we see that the functional I is bounded from below and the minimization 
problem 

r/ = inf {/(w)|w G W 1 ' 2 ^)} (4.22) 

is well-defined. 

Let {w^} = | fw^, . . . , Wm ^ | be a minimizing sequence of (|4.22p . It is easy to see that 
the function f(t) = ere* — rjt, where a,rj are positive constants, satisfies the property that f(t) — > 

00 as t — > ±oo. Thus, we conclude from (|4,2ip that |^.P|(i = l,...,m) are bounded. As a 

result, j (i = 1, ... ,m) are bounded. Therefore, the sequences (i = 1, ... , m) admit 

convergent subsequences, which are still denoted by (i = 1, ... , m) for convenience. Then, 

there exist m real numbers Wi, ■ ■ ■ , G R such that wf — > wn (i = 1, . . . , m), as £ — > oo. 
Using (|4.21|) again, we infer that |Vio^ 1 (i = 1, ...,m) are bounded in L 2 (f2). Therefore, 

it follows from the Poincare inequality that the sequences j (i = l,...,m) are bounded 

in W l,2 {Q). Consequently, the sequences = l,...,m) admit weakly convergent sub- 

sequences, which are still denoted by = l,...,m) for convenience. Then, there ex- 

ist m functions u^°°' 1 £ W 1,2 (Q,)(i = l,...,m) such that w!f — > w^ 00 ^ weakly in W l,2 {£l) as 

1 ->■ oo (i = 1, . . . , m). Of course, ui- 00 -* G PF 1,2 (0) (i = 1, . . . , m). 

Set w^ 00 "* = 2^°°^ + ^j 00 '' (i = 1, . . . ,iti), which are all in W 1,2 (£l) naturally. Then, the above 
convergence result implies xdf — > w^ 00 ^ (i = 1, . . . , m) weakly in TU 1,2 (il) as t — > oo. Noting that 
the functional /(w) is weakly lower semi-continuous in W 1,2 (£l), we conclude that (u^ 00 "*, . . . , w^^) 
is a solution of the minimization problem (I4.22p and is a critical point of /(w). As a critical point 
of I(w), it satisfies the equations fl4. 131) — (14. 15f) . 

Since the matrix R is positive definite, it is easy to check that /(w) is strictly convex over 
W 1 ' 2 ^!). As a result, the functional /(w) has at most one critical point in W 1,2 (Q), which implies 
the uniqueness of the solution to the equations ()4.13p - (j4.15p . 
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As in |35| . we briefly remark that we can also find a critical point of the functional 7 by a 
constrained minimization procedure. 

To proceed, we rewrite the constraints (|4.8p as 



Ji(w) = / e< +v * dx = Ki, i = l,...,m. (4.23) 

n 



By the assumption (|2.63p . we have Ki > 0,i = l,...,m. Then we consider the constrained 
minimization problem 

T] = inf {/(w)| w G W 1 ' 2 ^) and satisfies flOg]) }. (4.24) 

If the problem (|4.24p has a solution, say, w = (wi, . . . ,w m ), then there exist some numbers (the 
Lagrangian multipliers) fj,%, . . . , [A m S M such that 

(D(I + f i 1 J 1 + fi 2 J 2 + --- + fi m J m )(^))( W ) = 0, VweF 1 ' 2 ^). (4.25) 

Using a series of test configurations w 1 = (5n, . . . , Si m ),i = 1, . . . , m in (|4.25p successively, we have 

Ln/ii + L 21 fi 2 = 0, ) 
Lufii + L l+ ufi i+1 = 0, i = 2, . . . , m - 1, > (4.26) 

Lmml^m = 0, J 

which imply Hi = [i 2 = ■ ■ ■ = )J- m = 0. This is to say that, the constraints do not lead to the 
undesired Lagrangian multiplier problem, and any solution of the constrained minimization problem 
(|4,25p is a critical point of the functional ()4.18|) itself. 
From the constraint (14.23p . we see that 

Ul = Ln^ = \nK 1 - In / e u ° +Lll,bl dx, (4.27) 

Jn 

v { = Lu^w^ + LaWi = InKi - In / ^o+^-i^-i+^^dx, i = 2, . . . , m. (4.28) 

Jn 

From (I4.19p . we have 

J(w) - ttt E H V ^ll2 = E / ^^o+^dx - E / K^dx, (4.29) 
Plugging ()4.28p into (|4.29p . and using Jensen's inequality, we have 

m m m m 

iw - zx E ii V ^ ni ^ - E Ki ln Ki + E K * ln ai = E ^ ln ( 4 - 3 °) 

j=l i=l i=l i=l 

where <7j = |fi|exp^-p| J n u^dx^j ,i = l,...,m. From (|4.3Q[) we know that the functional I is 
bounded from below. 

Let {w^} = | \Wi\ . . . , Wm^j | be a minimizing sequence of the problem (j4.24j) . We conclude 
from (j43U1) that jwf 5 } (i = 1, . . . m) are bounded in W l,2 (Vt). Without loss of generality, we 
may assume |w^|(* = 1, ...m) converge weakly in W 1,2 (£l). The Trudinger-Moser inequality 
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(|4T2jl and (14271) and (14281) imply that {wfH (i = 1, . . . m) also converge. Therefore, the sequence 

| \ . . . , Wm^j | has a weak limit in W 1,2 (Q) as £ — > oo. Noting that the constraint functionals 
are weakly continuous and the functional / is weakly lower semi-continuous, we see that the weak 
limit of | . . . , Wm^j | must be a solution of the problem (|4.24|) . Since / is strictly convex, this 

solution is unique. Thus the constrained minimization procedure is carried out as well. 

Finally, integrating the equations (|42|) - (|44p or (|4.5|) over f2, we see that the quantized integrals 
(|2.64p are established. 

5 The special case where a = 

When a = 0, we denote the matrix R by R. Then 



/ 1 -1 
-13-2 
0-2 5 -3 



R 













-(m — 2) 2m — 3 — (m — 1 
-(m-1) 2m -1 J 



By a direct calculation, we see that the leading principal minors of R are 

Ri = i\, i = 1, . . . , m. 
Then by our formula (j3.6|) for the inverse of R, we have 



m ^ 

(R~ 1 ) ij = ^2-, j = i,i + 1, . . . ,m, i = l,...,m. 



that is 



ir 



/ m 


m 

El 


m 

E{ •• 


A 

m 


i=i 


Z=2 


Z=3 




m 

E} 


m 

E} 


m 

E{ •• 


l_ 

m 


1=2 


Z=2 


1=3 




m 

E} 


m 

El 


m 

E{ ■■ 


l_ 

m 




Z=3 


1=3 




\ m 


j_ 

rn 


J_ 

m 


. ±) 

m / 



From this expression of R 1 , we have 

m — i + 1, i = 1, 
Thus, applying Theorem 12.31 we obtain Theorem 12.41 



, m. 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 
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An interesting feature of the case a = is that the matrix-valued charge densities j° and J® 2 
may coincide to satisfy the relation [33] 

f = J°X2 = (5.6) 

so that, using (|2.48p . the minimum energy is seen to be directly related to the vortex numbers 
according to the quantization formula 



E = i fx \Q + 2R l2 \ = ^fx 



TV J j°dx + 2Tr / Jf 2 dx 

N-l 



kfJ* 

2^ 



Tr / B dx 



rO 



k/l ^2 (N ~ i) n i 



(5.7) 



i=l 



as stated in 1 33 1 . 



6 Some concrete examples 

When a > 0, the matrix computation quickly becomes rather involved for large N and the results 
are not as explicit as the case for a = 0. However, if iV is low and concrete, we can readily apply 
Theorem 12, 31 to obtain explicit results for the problem. As an illustration, we work out the equations 
for N = 3 (or m = 2) and N = 4 (or m = 3), respectively, as examples. 
We first consider the case when m = 2. The equations in this case are 

m 



Am = A([2a 2 + l]e Ul -[a 2 + l]e 112 -l)+4vr^5 PM , (6.1) 

s=l 

n.2 

Au 2 = X(-[a 2 + l]e Ul + [2a 2 +3]e M2 - l) +4vr^^ s , (6.2) 



3=1 

for which the associated coefficient matrix R reads 



2a 2 + 1 -(a 2 + 1) 
-(a 2 + l) 2a 2 + 3 



with the eigenvalues 



Ai, 2 = 2a 2 ± Va 4 + 2o 2 + 2 + 2. (6.4) 

Therefore 

A = 2 (2a 2 - \/a 4 + 2o 2 + 2 + 2) . (6.5) 
On the other hand, the inverse of R is 

«-'- 4 ', f 2 f +3 ° 2 2 +1 V (6.6) 

3a 4 + 6a 2 + 2 I a 2 + 1 2a 2 + 1 i y 1 



Thus, by definition, we have 



3a 2 + 4 3a 2 + 2 

n 3a 4 + 6a 2 + 2' r2 3a 4 + 6a 2 + 2" 1 j 



In view of Theorem 12.31 we obtain the following explicit results. 
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Theorem 6.1 For any a > 0, A > ; consider the equations (|6.ip - (|6.2p . 

(i) The equations (|6.ip - (|6.2p have a unique solution overM? satisfying the boundary condition 



3a 2 + 4 , 3a 2 + 2 

u\ — > in — k , U2 — > m ■ 



\x\ —> oo. 



3a 4 + 6a 2 + 2' 3a 4 + 6a 2 + 2 : 

Moreover, this solution obeys the following exponential decay estimate near infinity: 

2 

J2(ui(x) -Inn) 2 < C(e)e-^- e )^l x l, 



(6i 



(6.9) 



where e G (0, 1) is arbitrarily small, C(e) is a positive constant depending on e, r\,T2 are given by 
(|6.7p . and Ao is given (|6,5p . 

fiij T/ie equations (l6.ip -( [6T2j) over a doubly periodic domain £1 have a solution if and only if 
the following two inequalities 



4vr([2a 2 + 3]ni + [a 2 + l]n 2 ) < A|0|(3a 2 +4), 
4vr([a 2 + l]n 1 + [2a 2 + l]n 2 ) < A|0|(3a 2 + 2), 

hold simultaneously. Furthermore, if a solution exists, it must be unique. 

Next we consider the case when m = 3. The equations are 



3.10) 
3.11) 



ll 2 



A Ul = A ( [2a 2 + l]e ni - [a 2 + l]e U2 - l) + 4vr ^ 8 P1>S , 

8=1 

Au 2 = A (- [a 2 + l]e Ml + [2a 2 + 3]e" 2 - [a 2 + 2]e" 3 - l) + 4tt <5. 
Au 3 = A (- [a 2 + 2]e" 2 + [2a 2 + 5]e" 3 - l) + 4tt ^ A 



P2,s • 



J P3,s ) 



3.12) 
3.13) 
3.14) 



so that the associated coefficient matrix R takes the form 

/ 2a 2 + 1 -(a 2 + l) \ 



R 



(a 2 + 1) 



V 



2a 2 + 3 -(a 2 + 2) 
-(a 2 + 2) 2a 2 + 5 J 

We will not compute all the eigenvalues of R. But, instead, we write down the inverse of R: 

/3a 4 + 12a 2 + 11 2a 4 + 7a 2 + 5 a 4 + 3a 2 + 2 
2a 4 + 7a 2 + 5 4a 4 + 12a 2 + 5 2a 4 + 5a 2 + 2 
^ a 4 + 3a 2 + 2 2a 4 + 5a 2 + 2 3a 4 + 6a 2 + 2 



3.15) 



R- 



1 



4a 6 + 18a 4 + 22a 2 + 6 
which gives us 



(6.16) 



3a 4 + 11a 2 +9 
2a 6 + 9a 4 + 11a 2 +3' 

4a 4 + 12a 2 + 6 
2a 6 + 9a 4 + 11a 2 +3' 

3a 4 + 7a 2 + 3 
2a 6 + 9a 4 + 11a 2 +3' 
Applying Theorem 12.31 again, we can state 



>'3 



(6.17) 
(6.18) 
(6.19) 
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Theorem 6.2 For any a > 0, A > 0, consider the equations (|6.12p - (|6.14p . 

(i) The equations have a unique solution overM. 2 satisfying the boundary condition 

Ui — >• In rj, |x| — > oo, i = 1,2,3, (6.20) 

where ri,r2,r^ are defined by ()6. 17j) — (j6. 19p . Moreover, this solution satisfies the following expo- 
nential decay estimate at infinity: 

3 

^2(ui(x) - lnr-i) 2 < C(e)e-( 1 - £ )^ a; l, (6.21) 
i=i 

where e 6 (0, 1) is small, C{e) > depends on e, and Ao is twice the smallest eigenvalue of R. 

(ii) The equations ()6.12|) (j6.14|) over a doubly periodic domain Vt have a solution if and only if 
the following three inequalities 

2vr ([3a 4 + 12a 2 + ll]m + [2a 4 + 7a 2 + 5]n 2 + [a 4 + 3a 2 + 2]n 3 ) < A|ft|(3a 4 + 11a 2 + 9), (6.22) 
7r([2a 4 + 7a 2 + 5]ni + [4a 4 + 12a 2 + 5]n 2 + [2a 4 + 5a 2 + 2]n 3 ) < A|0|(2a 4 + 6a 2 + 3), (6.23) 
27r([la 4 + 3a 2 + 2]ni + [2a 4 + 5a 2 + 2]n 2 + [3a 4 + 6a 2 + 2]n 3 ) < A|0|(3a 4 + 7a 2 + 3), (6.24) 

hold simultaneously. Furthermore, if a solution exists, it must be unique. 

In general, for any N >2, since we have established the existence and uniqueness of a solution 
realizing n« arbitrarily prescribed zeros for each fi, i = 1, . . . , N — 1, we see that the solutions 
precisely depend on 

N-l 

2n = 2j2 n i ( 6 - 25 ) 

i=i 

continuous parameters which are the coordinates of the vortex points 

Pi,i, ■ ■ ■ ,Pl,m,- ■ ■ ,Pn-i,i, ■ ■ ■ ,PN-i,n N -i- (6-26) 

It may be relevant to recall the index-theory study of Weinberg [19] on the Abelian Higgs BPS 
equations which shows that the dimension of the moduli space of multiple vortex solutions is 2n. 
This number is precisely the number of the coordinates in M 2 which are needed to determine the 
locations of n vortices, in view of the classical existence and uniqueness results in |31j . There are sim- 
ilar studies in the context of non-relativistic Chern-Simons equations [33] and of superconducting 
cosmic strings [32]. It will be interesting to know the dimension of the moduli space of the multiple 
vortex solutions of the BPS vortex equations (|2.15p - (|2.17p of Kim-Kim-Kwon-Nakajima [33] in 
the ABJM theory. Thus our existence and uniqueness results indicate that, within the ansatz of 
Kim-Kim-Kwon-Nakajima |33j . the dimension of the moduli space of the n- vortex solutions of the 
BPS equations (|2.15p - (|2.17p is also 2n, although, generally when no specific ansatz is specified, it 
should depend on the gauge group index N as well [17] . 
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